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The relevance of intrinsic (or primordial) transverse momentum of partons in the inclusive pro- 
duction of particles at high energy and moderately large pr has been known for a long time, begin- 
ning with Drell-Yan and diphoton processes, and continuing with photon and meson production in 
hadronic collisions. In view of its renewed interest in the context of polarized processes and single 
spin asymmetries we perform, in the framework of perturbative QCD with the inclusion of spin and 
k± effects, a detailed analysis of several such processes in different kinematical situations. We show 
that the inclusion of these effects leads, at the level of accuracy reachable in this approach, to an 
overall satisfactory agreement between theoretical predictions and experimental unpolarized data, 
thus giving support to the study of spin effects and single spin asymmetries within the same scheme. 
We present results for transverse single spin asymmetries, generated by the so-called Sivers effect, 
in inclusive pion and photon production in proton-proton collisions. We compare our results with 
the available experimental data and with previous results obtained using simplified versions of this 
approach. 

PACS numbers: 12.38.Bx, 13.88.+e, 13.85.Ni, 13.85.Qk 



I. INTRODUCTION 



It has been known for a long time that the intrinsic (or primordial) transverse momentum of partons inside hadrons 
involved in high-energy processes may play a relevant role. A typical example is the Drell-Yan process, where 
the primordial, non-perturbative transverse momentum of the initial partons is directly related to the lower part 
of the transverse momentum spectrum of the observed lepton pair. A similar example is diphoton production in 
hadronic collisions. Since the first applications of perturbative quantum chromo dynamics (pQCD) to inclusive particle 
production in hadronic collisions, the role of these effects has been of some interest and several generalizations of 
the usual coUinear pQCD approach have been presented [1] (throughout this paper, by coUinear pQCD we mean 
the approach in which intrinsic transverse momentum, fcj^, effects are integrated out up to a factorization scale 
and all partons/hadrons are assumed to be coUinear with parent hadrons/partons). In recent years, several papers 
have reconsidered these effects, since coUinear pQCD, even at next-to-leading order (NLO) seems to underestimate 
experimental results for photon and pion production in hadronic collisions in the central rapidity and moderately 
large pt region [2-4]. It was found that the inclusion of intrinsic transverse momentum effects allows in most cases 
to reconcile theoretical calculations with experimental results. This requires a relatively large average transverse 
momentum, showing some dependence on the cm. energy of the process considered. This probably indicates that 
some effects due to higher order pQCD corrections are effectively embodied in the intrinsic momentum contributions, 
a point which needs to be further clarified. 

Almost independently of the above mentioned studies, the role of intrinsic transverse momentum has received a lot 
of attention in the context of polarization effects in inclusive particle production at high energy and moderately large 
Pt, in particular regarding transverse single spin asymmetries (SSA). Contrary to former expectations of pQCD [5], 
there are several experimental observations showing that SSA can be large in this kinematical regime. It was originally 
suggested by Sivers [6] that pQCD with the inclusion of transverse momentum effects in parton distribution functions 
could be able to explain these results. Sivers' suggestion has been further extended in a number of subsequent papers 
[7, 8]. There are several alternative possible explanations for SSA, basically non-perturbative in nature. Recently, 
Bourrely and Soffer [9] have claimed that most of the experimental data on SSA cannot be explained by pQCD, on 
the basis that the coUinear approach fails to reproduce the corresponding unpolarized cross sections (which enter 
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the denominator of the SSA) by 1-2 orders of magnitude. As a matter of fact, in all previous papers on SSA using 
pQCD-based approaches, a detailed analysis of unpolarized cross sections was not performed. The main reason for 
this was that a number of effects (scale dependences, NLO corrections) might cancel out in the SSA, which are ratios 
of (polarized) cross sections. Moreover, data on unpolarized cross sections in the kinematical regions where SSA are 
measured to be relevant are scarce. Considering the increasing theoretical and experimental interest in high-energy 
spin effects, a detailed combined analysis of unpolarized cross sections and SSA, within the same pQCD approach, is 
timely and worthwhile. 

Based on these considerations, in this paper we first perform a detailed analysis of a large amount of data on 
unpolarized cross sections for inclusive particle production in several processes and in different kinematical situations. 
Let us stress that our approach and aim are different in many respects from (and complementary to) those in Ref.s 
[2-4] : our main interest remains the study of SSA. Most of the interesting results on SSA are in kinematical regions 
different from those preferentially considered in the above mentioned papers, that is at 1 GeV/c < Pt < 3 GeV/c 
and large \xf\ = 2\pl\/^/s. Moreover, calculations of SSA are at present limited to the inclusion of leading-order 
contributions. Therefore, in this paper we limit ourself to study unpolarized cross sections at the same level of 
accuracy. In fact, our goal is not that of reproducing with very high precision the experimental results on unpolarized 
cross sections. Rather, wc want to show that the same approach used for SSA, at the same level of accuracy, is able 
to reproduce unpolarized cross sections for several processes and different kinematical situations, up to a reasonable 
factor (of the order of 2-3) which may be attributed to NLO corrections (the so-called iiT-factors), scale dependences, 
etc., which we expect, at least, to partly cancel out in the SSA. 

To this aim, we use here a generalization of the usual coUinear pQCD approach to inclusive particle production with 
the inclusion of spin and transverse momentum effects. We will take into account k± effects for all partons involved in 
the process, including the fragmentation process where required, which has been neglected in most of the Ref.s [2-4]. 
This considerably improves the treatment of k± effects with respect to Ref.s [7, 10], where a simplified approach was 
adopted, including k± effects only at first non-vanishing order and with a simplified partonic kinematics. Therefore, 
in this paper we will reconsider SSA for several processes already analyzed in Ref. [7], comparing the results obtained 
in the two approaches. 

The plan of the paper is the following: in section 11 we summarize our approach, which will be largely utilized in the 
rest of the paper. In section III we present and discuss our results concerning the unpolarized cross sections for the 
Drell-Yan process and for inclusive pion and photon production in hadronic collisions. A detailed comparison with 
available experimental results in several kinematical configurations is given. In section IV we then present our results 
for SSA in inclusive pion and photon production in hadronic collisions, including only the so-called Sivers effect [6-8] 
as the possible mechanism for generating the SSA. Our conclusions are presented in section V; the Appendix collects 
useful relations regarding the kinematics. 



II. PHENOMENOLOGICAL APPROACH 



In this section we present the basic ideas, ingredients and formulae of our approach; details on kinematics are given 
in the Appendix. Let us consider the inclusive production of large px hadrons in high-energy hadron collisions, that 
is the process AB C X. The starting point of our approach is the well-known collinear pQCD factorized expression 
for the corresponding differential cross section: 

^3^^ = XI J dXadXbdz fa/A{Xa,Q'^) fb/B{Xb,Q'^) ^| 6{s + t + u) Dc/c{z, Q"^) ■ (1) 

Here Xa, Xb (z) are the light-cone momentum fractions of partons a, b (hadron C) with respect to the parent hadrons 
A, B (parton c); /(x, Q^), D{z,Q^) are the well-known non-perturbative parton distribution functions (PDF) and 
fragmentation functions (FF) respectively, where the intrinsic transverse momentum dependence has been integrated 
over up to the factorization scale Q^. Furthermore, s, t, u are the Mandelstam variables for the partonic scattering 
process, ah cd, with all (massless) partons taken to be collinear with the corresponding hadrons. The Dirac delta 
function accounts for energy-momentum conservation in the elementary process and allows to express e.g. z as a 
function of Xa and Xb and the hadronic Mandelstam invariants (see also the Appendix). 

Extending this formalism with the inclusion of intrinsic transverse momentum of partons is not trivial at all and 
poses several problems. To start with, a complete formal proof of the factorization theorem is still missing in this 
case. The validity of factorization has been conjectured by Collins [11] and proved only for the Drell-Yan process and 
in e+e~ annihilation [12]. Quite recently a factorization formula for semi-inclusive deep inelastic scattering in the 
current fragmentation region has been proved [13]. A formally complete definition of transverse- momentum dependent 
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PDF's and FF's and of their evolution properties is still missing. In a partonic approach and at leading twist in the 
factorization scale Q the usual coUinear parton densities fa/Ai^a^ Q^) ^re simply generalized to transverse- momentum 

dependent distributions /a/yi(a;a, fe-Lo, Q^), where A;_lo is the transverse momentum of parton a with respect to the 
light-cone direction of hadron A, such that 

fa/AiXa.Q^) = j<fk±afa/A{Xa,k±a,Q'^). (2) 

Analogously, the usual parton fragmentation fimction Dc/d^, Q^) generalizes to Dc/,,(z, k±C; Q'^), where k±c is the 
transverse momentum of the final hadron C with respect to the light-cone direction of the fragmenting parton c, and 

Dc/c{z, Q^) = J d'^k^c Dc/c{z, k±c, Q') • (3) 

Intrinsic transverse momentum effects are of higher twist; consistently, a complete higher- twist treatment of the 
process would be required. However, wc arc far from being able to perform this complete analysis, since it intro- 
duces new, unknown non-perturbative PDF's and FF's and quark-gluon correlations; furthermore, the usual partonic 
interpretation of PDF's and FF's would be lost in this case. 

Despite these problems, there is a strong phenomenological evidence that transverse momentum effects in the 
partonic kinematics and dynamics are very relevant and may contribute to account for several experimental results 
for (un)polarized cross sections in inclusive particle production, difficult to explain in the collinear pQCD approach. 

This is why this generalized approach has been extensively employed in recent years; some of the resulting phe- 
nomenological outcomes are quite promising, although much work remains to be done on its formal aspects. This 
leads to a modified expression for the diff'erential cross section, Eq. (1), which now reads 



Ec da 



AB^CX 



d^pc 



^ / dXad^k±a dXbcfk±b dzcfkc S{kc ■ Pc) fa/A{Xa, k±a, Q^) fb/B{Xb, k±b, Q^) 



a,b,c.,d 



X 



{xa,Xb,s,t,u)-8{s^t^u)^J(z,\kc\)Ocic{.z.,kc-,Q )• (4) 



XaXbS dt TT Z 



Notice that the partonic cross sections da/dt depend, via the Mandelstam variables, on the intrinsic transverse 
momenta of partons. Notice also that the integration over the transverse momentum of the observed hadron C with 
respect to the light-cone direction of the fragmenting parton c runs over a generic three- momentum kc, defined in 
the cm. of the initial hadrons A and B (we will call this the hadronic cm. frame in the following); the delta function 
S{kc ■ Pc) guarantees that kc is in fact always orthogonal to Pc, the parton c three-momentum, and also the correct 
normalization of the fragmentation function, according to Eq. (3), which makes kc completely equivalent to k±c- 
The extra phase-space factor J(z, \kc\)/z'^ is the proper invariant Jacobian factor for the transformation from the 
parton momentum p^ to the hadron momentum p^ with inclusion of transverse momentum effects; the term s/{xaXbs) 
restores proper flux factor for non-collinear (in the hadronic cm. frame) parton-parton collision (see Appendix). 

Let us stress at this point that, once collinear pQCD is complemented with the inclusion of fc^ effects, not only 
Eq. (1) modifies to Eq. (4), but additional contributions to the unpolarized cross section are in principle possible 
[14-16]. These contributions to the unpolarized cross section could play a relevant role when considering asymmetries 
with respect to some measurable azimuthal angle related to the partonic process (e.g. in the Drell-Yan process, see 
[15]); however, we have explicitly checked that they give negligible contributions to the cross section for the inclusive 
process AB C X, where all "internal", partonic variables arc integrated over. Therefore, they are of little relevance 
for this analysis, and will not be further considered in the rest of the paper. A more complete discussion will be given 
elsewhere [16]. 

Eq. (4) has been widely used in the literature, starting from the pioneering work of Fcynman, Field and Fox [1]. 
Several papers have shown in the recent past that intrinsic fc^'s are indeed necessary in order to improve the agree- 
ment between next-to-leading order (NLO) pQCD calculations and experimental data for inclusive pion and photon 
production at high energies and moderately large px [2-4] . 

An expression similar to Eq. (4) holds also for the numerator of a SSA, {da^ — da^)/{da^ + da^) cx dA^a/da, 
replacing, for the polarized particles involved, the corresponding unpolarized PDF, FF and partonic cross sections by 
the appropriate polarized ones, A^f, A^D, dA^a (see Ref. [10] and section III for more details). At leading twist 
there are four such new spin and k± dependent functions to take into account: 

^'^fa/A-' = fa/A-! {x , k±) - f^/Al {x , k±) ; A^/aT/A = fal/Aix, k±) - fal/A{x, k±) ] (5) 

A^^c/cT = Dc/ct{z,k±) - Dc/ci{z,k±) ; A^L»ct/c = Dct /c{z,k±) - Dci/c{z,k±) , (6) 
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two in the PDF sector, Eq. (5), and two in the FF sector, Eq. (6); the first functions in Eq.s (5), (6) arc respectively 
the so-called Sivers [6] and Collins [11] function. The second ones are respectively the function introduced by Boer 
and Mulders [14, 15] and the so-called "polarizing" FF [14, 17]. For a general overview on spin and fe_L-dependent 
PDF's and FF's see also Ref. [18]. 

Regarding SSA, in this paper we will focus on the Sivers effect in inclusive pion and photon production. For a 
phenomenological study of the role of the Collins effect within the same approach and at the same level of accuracy 
see [19]. SSA and Sivers effect in the Drell-Yan process have been discussed in [20]. 

In the following sections, we will first concentrate on unpolarized cross sections, applying Eq. (4) to several processes 
in different kinematical situations. Some details of the approach may differ for different processes and will therefore 
be discussed in the appropriate sections. We now briefly comment on how some more general ingredients required in 
practical calculations based on Eq. (4) have been fixed. 

1) Concerning the fc^-dependent PDF's (FF's) we assume that the x (z) and k± dependences factorize; we use 
a Gaussian-like, flavour-independent behaviour for the A;_L-dependent part; that is, neglecting for the moment the 
dependence on the factorization scale, we take {k± = \k±\) 

fa/A{x, k±a) = fa/A{x) ^ e'^^' ; Dc/c{z, k±c) = Dc/c{z) — e-^'" '^ic , (7) 

TT TT 

where the parameter /? (/?') is related to the average partonic (hadronic) k±_ by the simple relation l//3(/3') = 

lh.2 \l/2 
\'^JLa(C)l ■ 

2) The cross sections for the elementary, partonic scattering, da"-^^"'^ /di, are calculated at leading order in the 
strong coupling constant power expansion, including in the partonic kinematics the full dependence on the intrinsic 
transverse momenta k±. That is, the partonic cross sections in Eq. (4) will depend on the properly fc^-modified 
partonic invariants s, t and u. At relatively low px the inclusion of k± dependence might make one or more of the 
partonic Mandelstam variables become smaller than a typical hadronic scale. This configuration would correspond 
to a situation where the propagator of the exchanged particle in the partonic scattering becomes soft; the observed 
transverse momentum is thus generated mainly by fluctuations in the intrinsic fc_L distribution and not by the hard 
scattering. In this case perturbation theory would break down. In order to avoid such a problem and extend this 
approach down to px around 1-2 GeV/c (where most of data on single spin asymmetries have been collected) different 
ways have been proposed. Following [2], we introduce a regulator mass, fi = 0.8 GeV, shifting all partonic Mandelstam 
variables, that is we take 

i^i—jjL^, u^u — jjL^, s^s + 2ij?. (8) 

3) Another somehow related source of potential ambiguity is the behaviour of the strong coupling constant, as{Q^), 
in the low regime. We adopt the prescription originally proposed by Shirkov and Solovtsov [21], using for as the 
expression 



Po 



1 A2 



log(0VA2) A2-Q2 



(9) 



where as usual /3o = (33 — 2n/)/127r, n/ being the number of active flavours (we use n/ = 4), and A = 0.2 GeV/c. 
According to Eq. (9), at large q,^ reduces to the standard LO expression, while at low its behaviour is well 
under control without the introduction of any extra parameter, like a freezing scale parameter or a dynamical gluon 
mass. We have explicitly checked that other prescriptions give similar results. 

Note that for the Drell-Yan process Eq. (4) simplifies a lot (see section II A); moreover, the large scale involved 
(namely the invariant mass of the observed lepton pair) removes the danger of critical kinematical regions discussed 
above. 

4) Another important ingredient is the choice of the factorization scale Q (we use one single scale for renormalization 
and factorization scales) governing the pQCD evolution of PDF's and FF's and entering the QCD running coupling 
constant, a.^- It is well known that LO calculations, and in some cases NLO too, can strongly depend on this choice. 
Moreover, there is not in the literature a unique prescription for the choice of the factorization sc;ale Q in inclusive 
hadron production, and several possible alternatives have been suggested. Typical scales adopted are the observed 
hadron transverse momentum, pT^ or pr/'^, and the transverse momentum of the fragmenting parton in the partonic 
cm frame, p^p, or p^/2. Apart from the Drell-Yan process, where the natural scale is the invariant mass of the 
produced lepton pair, throughout this paper we will adopt the factorization scale Q = p^/2. 

5) Theoretical results for cross sections at LO and NLO are usually compared to get an estimate of the uncertainty 

related to the choice of the factorization scale and indications on the convergence of the perturbative expansion. 
The ratio of cross sections evaluated at NLO and LO respectively gives the so-called if-factor for the process under 
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consideration. It is often assumed [2] that this ratio is independent of intrinsic transverse momentum effects, so 
that it can be estimated in the usual coUinear pQCD approach. Throughout this paper we will conform to this 
prescription; that is, in comparing our LO results for cross sections with the corresponding experimental results, 
we will include coUinear pQCD iiT- factors evaluated by means of independent numerical codes [22-24]. For the 
processes under consideration in this paper, the i^-factors show some general behaviour with respect to the relevant 
kinematical variables [24] : they decrease smoothly with the increasing of the observed hadron transverse momentum, 
Pt (the variation being steeper at the lower pt considered), at fixed ^/s; they also decrease mildly with -^i at fixed 
Pt, whereas they show a tiny dependence on x^. For simplicity and clarity, in our plots of cross sections we will 
always show our numerical curves rescalcd by a fixed (for that curve) i^-factor, obeying the general trends of above 
and chosen according to the results of the numerical code INCNLL of Ref.s [22, 23]. We think this should help in 
clarifying the role played by ii"- factors in our results. The mild variations of the computed ii'-factor for each point of 
a given curve should not change our general conclusions. More details on the estimate of /C-factors are given in the 
following sub-sections. 

Let us finally recall that our main aim is not to perform an overall best fit to the available data (to this end, a 
full NLO approach would be required). Rather, we want to show to what extent the same approach adopted in the 
description of SSA is able, when complemented by coUinear pQCD i^-factors, to give a reasonably good account of 
a large set of cross-section data for different processes in different kinematical situations. We feel that the approach 
here described is reliable in this respect. Notice also that, as we have explicitly checked, the main interest of our 
study, namely single spin asymmetries, are not strongly affected by the choices discussed above. 

In the following sub-sections we discuss separately and in more detail the processes under study. 



A. Drell-Yan process 

Let us now consider lepton pair production in pp, pp collisions, the Drell-Yan process. In LO, coUinear pQCD the 

lepton pair has no transverse momentum, q^i, in the hadronic center of mass frame. On the other hand plenty of 
data show a well-established qr spectrum, with an exponential behaviour at low qT (^2 — 3 GeV/c) turning into a 
power-like one at larger transverse momentum. NLO corrections and soft gluon radiation are important contributions 
both for the normalization of the cross section and in order to reproduce the observed power-like behaviour at large 
qr- However, to describe the behaviour at very small qr one has to include the intrinsic parton momentum, k±. 

Therefore, within our approach the Drell-Yan process can be considered a useful tool to gain information on the 
intrinsic transverse momentum of partons inside the initial colliding hadrons. 

Analogously to Eq. (4), the cross section for the production of an ^+^~ pair in the collision of two hadrons A and 
B (there is no need for any fragmentation function here) reads: 

dcr = d^k^a dXb d'^kj_b fa/A{Xa,k±a) fb/B{xb, kj_b) — - — da"''^'^^ , (10) 

J XaXbS 

ab 

where the elementary cross section da for the process a{pa) b{pb) — > £'^{p+) i~{p-) is given by: 

The differential cross section da depends on the variables 

s^M^ = {p,+Pbf^q\ y=iln^O±^, q^, (12) 

that is the squared invariant mass, the rapidity and the transverse momentum of the lepton pair; qo, qj, and qL 
are respectively the energy, transverse and longitudinal components, in the A-B cm. frame, of the four-vector 
q = Pa + Pb — P+ + P-- Using the relations: 

^=d^p-5{pl) p-=q-p+ dM^dy = 2dqodqL, (13) 

Eq. (10) can be written as 

^ ^ =^ [dXad'^k_LadXbd'^k±bfa/A{Xa,k±a) fb/B{Xb,k±b)S'^iPa+Pb- q) —^^f , (14) 



dy dM^ d^qr 

ab 
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where a^^ is the total cross section for the ah ■ 



'Lab 
On 



1 



2E+ {2tt)'^2M'^ 



£'^£ process: 
1 



^{{q - P+) ) I Mab^c+c- {p+,q) I 



(15) 



For the kinematical regimes of interest in this paper the dominating elementary contribution to the Drell-Yan 
process is the lowest order electromagnetic interaction, ^ 7* — > £'^£~, so that a,b = q,q with q = u,u, d, d, s, s and 



.qq _ 47ra^ 
~ 9M2 

The Dirac-5 function accounting for energy-momentum conservation in Eq. (14) contains the factors 

1 
2 
1 



(16) 



■ S{Ea + Eb 



{Xa + Xb)^ + 



qo) S{pza +Pzb 
2 



_|_ 

XaS X})S 



{Xa-Xb)—- 



A-2 

XaS 



A-2 

XbS 



In the following we shall only consider kinematical regions such that: 



4 « kl^^b 
where Eq. (17) simplifies into the usual collinear condition: 



■ 6{Ea + Eb- qo)6{pza + Pzb - qh) 



M 



Xb 



M 



(17) 



(18) 



(19) 



The Gaussian shape adopted for the k± dependent PDFs', Eq. (7), together with Eq. (19), allow us to analytically 
perform the integrations in Eq. (14) which becomes 



00 P'P' 



dydH'^d^qj. tts /32 + /32 



exp 



5Z ^9 fl/p^^o) fq/pi^b) ) 



(20) 



where P, (3 refer to the quark(antiquark) PDF respectively. By direct comparison with data we can then extract the 
/3, P parameters. The numerical values obtained in this way are of course related to the set of x dependent PDF's 
adopted; throughout this paper wc use the MRSTOl set [25]. In principle /? and /? could be x and flavour dependent, 
but in first approximation a reasonable description of data can be obtained by neglecting these dependences. In this 
case, ^ = /3 and l/p'^ = {q^)/2. 

There is a clear experimental evidence that the average transverse momentum (qr) of the lepton pair increases with 
the cm. energy, in agreement with pQCD calculations [26]. In our LO approach, see Eq. (20), this behaviour can be 
obtained by using an effective value of P decreasing with cm. energy. This way, one should more correctly interpret 
the parameter P as representing the effects of the primordial transverse motion of partons plus a component coming 
from NLO corrections. Since our first aim here is to extract information on the primordial intrinsic momentum, we 
write 



(Q't) — (9T)intr + (?T)pert j 



(21) 



where (gff')pert oc s, is a perturbative, energy dependent contribution, which will be estimated by comparison with 
data. A good description of data available in different kinematical regions can be obtained by choosing at the lowest 
cm. energies considered, ^/s ~ 20 GeV, 



( (g^)i„tr/2 )^/2 = l//3o - 0.8 GeV/c . 



(22) 



In Fig.s 1-5 we compare our estimates of the Drell-Yan invariant cross section in pp collisions, averaged over bins 
of the invariant mass M and at fixed rapidity, as a function of qr, with a collection of data from several experiments 
[27-29]. 

All curves are rescaled by proper if- factors varying in the range 1.5-1.8 in the different cases (see legend and caption 
of figures), in agreement with NLO calculations. The factorization scale is set to the lepton pair invariant mass. 

One can sec that experimental data are described very well for qx values up to 2-3 GeV/c, where (see Fig. 3) a 
power-like behaviour, entirely due to radiative effects, starts to set in. 
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B. Direct photon production 

Let us now consider inclusive prompt photon production in pp collisions (for a compilation of data and a complete 

and detailed NLO analysis in coUincar pQCD of prompt photon production in hadron-hadron interactions sec [30]). 

As pointed out in many papers (see for instance [31]) intrinsic transverse momentum effects can help in solving 
the discrepancy between experimental data and LO as well as NLO collinear pQCD calculations. In particular, what 
emerges [32] is that the steep px dependence of the measured differential cross section cannot be explained by any 
new improved PDF's (mostly for gluons); on the contrary the introduction of intrinsic transverse momentum effects 
significantly improves the comparison between theoretical calculations and experimental data. An attempt to give a 
more firm theoretical foundation, in terms of Sudakov form factors, of the phenomenological Gaussian fej_ smearing 
has also been carried out [33] . 

We employ here our LO pQCD approach (including proper /-T- factors from NLO collinear pQCD), neglecting 
possible photon fragmentation contributions, an issue related also to isolated photon cross sections measured at 
extreme high-energy collider experiments (not considered in our analysis) [30]. 

Starting again from Eq. (4), the invariant differential cross section for the process pp ^ jX then reads 



E^-^ — = (^^a(fk±adxbcfk±bfa/p{xa,k±a)fb/p{xb,k±b) S{s + i+u), (23) 

' ab 



where the basic partonic processes are the Compton process gq{q) — > jqiq) and the annihilation process qq — > 75. 

By exploiting the elastic constraint, s + t + u ^ 0, one of the integrations can be easily carried out, which fixes one 
of the light-cone momentum fractions (e.g. Xb); the 5-dimensional integral left is handled with the help of a VEGAS 
Monte Carlo routine which properly takes into account the full kinematics and all required kinematical cuts (see also 
the Appendix). 

The kinematical regime we are interested in is mainly the moderately large region of pr values (1-4 GeV/c), where 
a uniform kj^ smearing on a steep falling px distribution produces a significant enhancing factor. 

Similar studies (and conclusions) can be found in [3, 4] which however consider only photon production in the 
central rapidity region. Here we extend this approach to larger xp values, in view of its interest in studying SSA. 
A semi-analytical analysis with estimates of the enhancing factor resulting from the inclusion of intrinsic transverse 
momenta, in the full range of xf values, is under completion and will be published elsewhere [34]. 

Our main results are compared in Fig.s 6, 7 and 8 with a representative set of experimental data for pp collisions, 
both from fixed target experiments, E7G4 (FNAL) [35] and WA70 (CERN) [36], and from coUider experiments, R806 
(ISR) [37]. For fixed target experiments and cm. energies below 40 GeV we have adopted the same value of /3o 
as extracted in the Drell-Yan process, /3o = 1.25 (GeV/c)-i, corresponding to {kiy/"^ = l/(3o = 0.8 GeV/c. For 
larger energies (ISR), Fig. 8, a slightly smaller value of /?, /3 = 1 (GeV/c)^^ seems to better reproduce the low pt 
distribution. As said, we adopt the factorization scale Q — Py/2, where p^ is the photon transverse momentum in the 
partonic cm. frame. Notice that the non-collinear parton configuration in the hadronic cm. frame implies p^ 7^ pt- 

We use a NLO JT-factor which decreases mildly with px at fixed energy {K= 1.6 for 1 GeV/c < < 3 GeV/c 
and K = 1.3 for px > 4 GeV/c). This is consistent with direct calculations we performed using the numerical code 
by Aurenche et al. [22] for NLO order direct photon production. As an example, in Fig. 8, where for the sake of 
simplicity we adopted a fixed value, K = 1.5, use of a px dependent if- factor would imply an increase of the low px 
edge (where iv" ~ 1.9) and at the same time a reduction of the large px tail (where iiT ~ 1.2) in our estimates, leading 
to a slightly better agreement with the data than shown. 

For comparison in Fig.s 6, 7 and 8 we also show the corresponding LO, collinear pQCD results. As expected, at 
very large the intrinsic transverse momentum effects are negligible. In fact, fej, contributions behave like k_\_/px; 
moreover, at large px the spectrum is less steep than on the lower part, where even a small smearing produces a big 
effect. 

As it is clear from Fig. 7, collinear pQCD calculations may already give an accurate description of WA70 data 
{y/s = 23 GeV) . On the other hand, results including fej_ efi'ects compare with data equally well, while E704 data [35] , 
even if at a slightly lower energy {^/s = 19.4 GeV) and at lower px, are heavily underestimated without fe_L effects. 

Once more we think that the agreement between our calculations and experimental data is good enough for our 
purposes. 



C. Inclusive pion production 



Large px pion production in hadronic collisions is certainly the most interesting inclusive process from many points 
of view. In this case the discrepancies between collinear pQCD calculations and experimental data are even more 
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significant. This is especially true at cm. energies below 60 GeV, where theoretical calculations can underestimate 
experimental data by a factor 10 or larger, depending also on the choice of the factorization scale (see also [9]). Again, 
the inclusion of intrinsic motion of partons has been advocated in order to improve the agreement between theory 
and experiment. 

Most part of previous work has been addressed to inclusive pion production in the central rapidity region [2, 38]. 
On the other hand, as already mentioned, most of the interesting and puzzling experimental results on SSA are in 
the region of intermediate and large positive xp values; therefore a consistent study of unpolarized cross sections in 
this kinematical regime and within the same model is mandatory. 

To this end, we again turn to Eq. (4). The clastic constraint allows us to carry out explicitly the integration 
over z, see the Appendix for details. Again, we employ a VEGAS Monte Carlo routine to perform the remaining 
8-dimensional integral with proper kinematical cuts. 

With respect to the previous cases considered (the Drell-Yan process and prompt photon production) a new trans- 
verse momentum dependent function enters through the fragmentation process. Notice also that the unpolarized FF's 
are presently known with much less acciiracy than the nuclcon PDF's. In particular, all available parametcrizations 
(with the exception of the set by Kretzer, Leader and Christova (KLC) [39]) for the pion FF's are based on e+e~ 
data for charged pion production, which do not allow to explicitly separate between 7r+ and tt" case; this separation 
can only be made under further assumptions, which remain to be tested. The z and A;_l dependences in the FF are 
chosen according to Eq. (7). 

Concerning the factorization scale, Q, as for the case of prompt photon production we use Q = where 
is the transverse momentum of parton c in the partonic cm. frame, that is = [ui/sY^^ (notice that in prompt 
photon production the outgoing elementary particle at LO is the photon itself). 

1. Neutral pion production 

Let us start considering inclusive neutral pion production. Using the fcx-dcpendcnt PDF's as inferred from the 
analysis of the Drell-Yan process and comparing our results with a large selection of cross-section data we can extract 
the (flavour independent) /3' parameter entering the pion FF, see Eq. (7). In this sub-section we will indicatively 
use the set by Kniehl, Kramer, and Potter (KKP) [40]. An optimized choice of /3' seems to favour an explicit z 
dependence, which somehow spoils the simple factorizcd form of Eq. (7), 

1 / P' (z) = { kl^{z) y =^0.7 z^'^il-zf-^^ GeY/c [KKP]. (24) 

Estimates of NLO corrections arc taken into account by rescaling our results by proper X-factors. With our choice 
of the factorization scale a iC-factor decreasing with energy and approaching unity at very large energies [24] is to 
be expected. This trend has been also checked by direct calculations performed with the help of the NLO numerical 
code for inclusive pion production by Aversa et al. [23]. 

In Fig. 9 we present our estimates of the invariant differential cross section for the process — > tt" X at different 
energies vs. pt and at central rapidity, y = 0. After rescaling each curve by the proper JC-factor, ranging from 2.5 to 
1, we get reasonable agreement with experimental data [41-44] over a large range of energies (from 20 to 200 GeV). 

A bit more problematic is the case of neutral pion production in the very forward region (very small scattering 
angles), as shown in Figs. 10 and 11. In particular for 6 < 10" — 15*^ even with the inclusion of intrinsic k± effects 
the discrepancy between theory and data at ISR energies [41], ^/s = 23.3 and 52.8 GeV, is still large. Nonetheless in 
these plots the LO collinear calculations (rescaled by the same JsT-factors) show how important the enhancing factor 
due to k± contributions can be. 

In Fig. 12 we compare our estimates, rescaled by a isT-factor K = 2.5, with data at ^/s = 19.4 GeV [42]; here fair 
agreement can be reached over a broad region of xp away from 0, and at different pt values. To our knowledge, these 
data are the only ones available covering almost the same kinematical region of the well-known E704 results on pion 
SSA [45]. 

Recent data from STAR [46] in comparable xp and pr ranges but at much larger energies are shown in Fig. 13, 
together with our theoretical curves, rescaled by proper JC-factors. Let us remind that we use for each curve a fixed 
-R'-factor, whereas use of a p-j -dcpcndent if-factor would imply an increase of the low E^^ edge and a reduction of the 
large Ej: tail. In this case wc; overestimate the data, but still a reasonable description is obtained. Let us notice that, 
since data on pion SSA come mainly from E704 experiment, in performing our optimized choice for /3', we have given 
a slight preference to unpolarized cross-section data in the same kinematical region [42] . A smaller value of 1 //?' would 
lead to a better agreement with STAR data without spoiling the good description of FNAL experimental results [42] . 
Let us also recall, in this respect, that the results presented cover a huge spectrum of different kinematical situations 
at different cm. energies. 
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We have also checked that the pion FF set by Kretzer [47] with a corresponding proper choice of the /3' parameter, 

= (fcl^(^))V2 = 1.421-2 (1-^)0.35 Qgv/c [Kretzer], (25) 

gives similar results. For the processes under consideration here the KLC set [39] is practically equivalent to the set 
ofRef. [47]. 

2. Charged pion production 

The study of inclusive charged pion production needs some additional caution; as it was already mentioned the 
available pion FF sets come from fits to experimental data on hadron production in e'^e^ annihilation, where only the 
sum of charged pions is detected. In order to get flavour (or charge) separation one therefore has to rely on some extra 
assumptions: for instance Kretzer [47] imposes, besides isospin symmetry, the following behaviour of the non-leading 
quark FF's at the input scale (/Uq)- 

Dn+/diz,fil) = {1- z)D^+/^{z,nl), (26) 
D„+/s{z,fio) = D„+/d{z,fJ-l)- (27) 

Under proper pQCD evolution, a FF set for tt" and tt^ can then be extracted. In other paramctcrizations, like the 
KKP set [40], no additional flavour separation assumption is made and only a FF set for neutral pions is available. 

Clearly the relative weight between valence (or leading) and sea (or non-leading) quark FF's can play a relevant 
role in determining the relative yields of 7r+ and Tr~ . 

For these reasons, we model two charged pion FF sets, starting from the neutral pion KKP parameterization. 
In one case (KKP-1) we assume that the non-leading (sea) quark FF's into a charged pion are all equal to the 
strange quark FF into a neutral pion: this leads to a strongly valence-dominated fragmentation mechanism at large z 
values. A second case (KKP-2) is defined along the lines of Kretzer set. More precisely, assuming isospin symmetry, 
-^71— /« = ^7r+/d and charge conjugation invariance, D„-/y_ = D„+ /^i, we define for KKP-1 

D^yaiz) = D„o/,{z) , (28) 

that implies 

Dn+/u{z) = 2D^o/„(z) - D^o/s{z) 
For KKP-2 set instead we impose 

D^^/diz) = (1 - zyD^y^iz) , (30) 

that implies 

D^+/^{z) = 2D„o/^{z)/[l + il-zr] (KKP-2), (31) 

where the parameter o, chosen reasonably around 1-1.5, will be fixed to a = 1.3 in the following. Clearly, these sets 
leave unchanged, and equal to those of the KKP set, the FF's for neutral pions. 

Let us remark that the analysis of Ref. [39], based on a simultaneous fit to e+e~ and semi-inclusive DIS results 
which in principle allows flavour separation, indicates a relatively small suppression of non-leading vs. leading quark 
FF's, as already assumed in Kretzer set. 

In Fig. 14 we show our estimates for the ratio of tt^ and 7r+ yields and compare them with FNAL experimental 
data [48] for three different energies and at scattering angle 6 = 77 mrad (at these energies, this corresponds almost to 
xf = 0) as a function of pr- Results using the three different options for pion FF's discussed above (Kretzer, KKP-1 
and KKP-2) are shown. Since all pion FF sets adopted assume that the FF for a neutral pion coincides with the 
average of those for charged pions, our estimates for tt" cross sections, Fig.s 9-13, together with the results of Fig. 14, 
give complete information for the charged pion sector. Therefore, we do not explicitly show our results for charged 
pion cross sections, which are in good agreement with data. 

From the ratio of 7r~, 7r+ yields wc sec that both Kretzer and KKP-2 sets give a reasonably good description of 
the experimental data, while KKP-1 set seems to underestimate the data, in particular at large px- Let us however 
remark that the KKP-1 and KKP-2 sets are model sets based on some assumptions which remain to be tested. We 
propose them here as representative of two very different behaviours with respect to the relative weight of leading and 
non-leading pion FF's. As such, our results are not intended to give preference to any set of FF's over the others. 



(KKP-1) . 



(29) 



10 



III. TRANSVERSE SINGLE SPIN ASYMMETRIES AND SIVERS EFFECT 



We come now to the other main subject of this paper. Parton intrinsic transverse momentum effects, as we have 
seen in the previous sections, can play a crucial role in reducing the discrepancies between pQCD estimates and 
experimental data for unpolarized cross sections in inclusive hadron production. What is more important from our 
point of view, they can be an essential ingredient in understanding some spin phenomena in the framework of pQCD. 
The large SSA observed in many reactions like A^B^CXotAB^ X (where H is typically a hyperon) can 
be described by certain spin-dependent effects generated by soft mechanisms in the presence of intrinsic transverse 
momentum. These effects and several applications to different processes, including various possible mechanisms of 
this sort, have been discussed in a series of papers [7, 10], keeping only dominant k± contributions and under some 
simplifying assumptions about parton kinematics and fej. distributions. 

Let us then consider the transverse single spin asymmetry 

Aj^ = '^'^^ ~ '^'^^ (32) 

where da^'^ stands for the invariant differential cross section for the process A^'^ + B ^ C + X , and t, J, indicates the 
transverse direction with respect to the production plane (corresponding to the ±Y direction in a reference system 
where the initial polarized hadron moves along the +Z direction and the observed hadron lies in the +X-Z half-plane). 
As it was already mentioned, at leading twist there are three possible soft mechanisms (and three corresponding new 
spin and k± dependent PDF's/FF's) contributing to this process, see Eq.s (5), (6), often referred to as "odd under 
naive time reversal" : 

a) The Sivers effect (and distribution function) [6], corresponding to the possible azimuthal dependence (around the 
light-cone direction of the parent nucleon) of the number density of unpolarized quarks (or gluons) inside a transversely 

polarized nucleon. 

b) The Collins effect (and fragmentation hmction) [11], corresponding to the azimuthal dependence (around the 
light-cone direction of the fragmenting parton) of the number density of unpolarized hadrons resulting from the 
fragmentation of a transversely polarized quark (or linearly polarized gluon). 

c) The Boer-Mulders effect (and distribution function) [14, 15], corresponding to the azimuthal dependence (around 
the light-cone direction of the parent nucleon) of the number density of transversely polarized quarks (linearly polarized 
gluons) inside an unpolarized nucleon. 

When considering the production of a transversely polarized, spin- 1/2 hadron (like e.g. in the process AB ^ X) 
another leading twist, T-odd effect has to been taken into account in the fragmentation process: 

d) The so-called polarizing fragmentation function [17], corresponding to the azimuthal dependence (around the light- 
cone direction of the fragmenting parton) of the number density of transversely polarized hadrons resulting from the 
fragmentation of an impolarized quark (or gluon). 

It should be noted that in the pure parton model, where partons are treated as physical, massless free particles, all 
these effects, when described as totally independent, factorizcid processes, vanish [49]. 

Let us remind that these functions may be defined also for gluons; in this case, however, instead of transverse quark 
polarization we must refer to linear gluon polarization [16, 50]. 

All of the above spin effects a) c) could contribute in principle to the process p ir -\- X , but, at the present 
stage, a complete phenomenological analysis is out of reach and perhaps not very significant, due to the relatively 
scarce experimental information available. On the other hand, it is very important to verify whether any of these 
possible contributions, when taken alone, may contribute and to what extent to the SSA observed, extracting even 
gross features of the corresponding PDF's/FF's. 

In the following we will then concentrate on SSA in inclusive pion and photon production generated by the Sivers 
effect [6] . SSA in Drell-Yan processes within the same approach have been considered in Ref. [20] . For a complementary 
phenomenological study along the same lines, but invoking only Collins effect, see Ref. [19]. 

A more formal and complete analysis including all possible leading-order, leading-twist contributions for single and 
double spin asymmetries in inclusive pion production is currently under way and will be presented elsewhere [16]. 

Let us then start from the most general expression for the number density of unpolarized quarks q (or analo- 
gously gluons), with light-cone momentum fraction x and transverse momentum fej_, inside a proton with transverse 
polarization P and three-momentum p. One can write 

/g/pT {x, k±) = fq/p{x, k±) + ^ A^/g/pT {x, k±) P ■ {p X k±) , (33) 



which implies 



fq/p{x,k±) = hfq/pT{x,k±) + fq/pl{x,k±)], (34) 



11 



^^fq/p^ (2;, k±) = fg/pi {x, k±) - fg/pi {x, k±) = A^/g/pT {x , k±) sin(^p - (f>kj , (35) 

where (f>p and (pkj^ arc respectively the azimuthal angle of the proton (transverse) polarization vector and of 
fe_L. Eq. (34) gives the unintegrated, fcx-dependent unpolarized PDF, while Eq. (35) defines the Sivers PDF, 
A^/q/pT (a;, /c^). Another common notation for the Sivers function is f^{x,kj_) [14]; it should be noted that the 
two definitions are not completely equivalent, the exact relation being A^/^/pt {x, k±) = {—2k±/Mp) f^{x, k±). 

In order to perform numerical estimates and carry out a comparison with available experimental data we introduce 
a simple model for the Sivers function (see also Ref. [20] for more details), similar to what was previously done for 
the so-called polarizing fragmentation function in Ref. [17]. 

Analogously to the case of the unpolarized PDF's and FF's, Eq. (7), we consider a simple factorized form for the 
Sivers function: 

A^./,/pT {x, k^) = A^/,/pT {x) h{k^) . (36) 
In order to satisfy the positivity bound (see Eq.s (34), (35)) 

|A^/g/pT(a;,fci)| < 2fg/p{x,k^), \/x,k^, (37) 

we put 

A^/,/pT(x) = 2Mq{x) f^,p{x) (38) 
h{kx:) = H(A;x) — e-^'^^i , (39) 

and we simply need to require 

\Mq{x)H{ks_)\<l yx,k±. (40) 

We actually impose J\f and H to be separately smaller than unity in magnitude, by choosing simple functional forms 
and dividing each of them by its maximum value: 

M,{x)=N,x^^{l-xt^ ^ "^Jj , |iV,|<l (41) 

W(fci) = V2e(a2-/32)fci exp [-(a^ - /S^) fci ] , a > (3 . (42) 
Notice that 'H{k±) must vanish at least like k± for k± — > 0. Eq.s (39) and (42) imply: 

A^/,/pT(a;,fci) = 2Mg{x)fg/p{x) ^ v'2e (a^ - /32) fc^ e-"''=i . (43) 

Since the positivity bound requires a > P, one can write = where r is a positive parameter smaller than 

one. We will keep this last parameter flavour and x-independent. 

Let us stress again that in principle one can also define a gluon Sivers function, for which the same expressions 
would be valid. In order to reduce the number of free parameters and assuming in first approximation that transverse 
spin-fej_ effects are valence-dominated (certainly a reasonable assumption for particles produced at large positive 
values of xf), we will restrict here to consider valence quarks, neglecting (for what concerns the Sivers effect only) 
possible contributions from sea quarks and gluons. We have explicitly checked that, for the processes under study, 
even maximizing the corresponding Sivers functions these contributions could play a role only at low xp, where 
is almost negligible (see also Ref. [51]). 

Under these assumptions we are left with seven free parameters which completely fix the Sivers function as given 
in Eq. (43); six for Afq{x) in Eq. (41), with q = u,d, plus r = in Eq. (42). 

A. Sivers effect in p^p — > nX 

By comparing our calculations for SSA in pion production with available experimental data we now try to fix the 
set of parameters entering the Sivers function. 
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According to Eq. (4), the numerator of the SSA, Eq. (32), in terms of Sivers effect alone reads 
Ett^ E^-^— = X] Jdxad'^k±adxbd'^k±bdz(fK6{K-p^)A^fa^/pr{xa,k±a]Q^)fb/p{xb,k±b;Q^) 

^ a,b,c,d 

s da"^^"'^ , , ' s . .,1 



{xa,Xb,s,i,u)-5{s + i + u)^J{z,\k.,,\)D^/c{z,kT,\Q'^). (44) 

The denominator of the SSA, da^ + da^ , is simply given by twice the unpolarized cross section, Eq. (4). 

Let us remind that we consider the kinematics of the reaction in the cm. of the colliding protons, where t,i 
indicates the transverse direction with respect to the production plane (corresponding to the ±y direction when the 
initial polarized proton moves along the +Z direction and the observed pion lies in the +X-Z half-plane). According 
to Eq. (35), in this frame = 7r/2, and 

fa/p^{Xa,k^a) = fa/pt {Xa, k±a) COs(0fe^ J . (45) 

As already discussed in section 11 C, energy-momentum conservation in Eq. (44) allows us to carry out e.g. the 
integration over z (see also the Appendix) leaving us with an 8-dimensional integral. This makes unpractical and 
very time-consuming a complete best-fit procedure over all seven free parameters. Moreover, a complete analysis 
of this kind would probably be premature and not very significant at this stage. Therefore, we here adopt a less 
rigorous approach and try to reproduce the experimental data, extracting an optimized choice for the Sivers function 
parameters from a limited set of choices in the available parameter phase-space. 

To this end, we consider the experimental data for pion SSA of the E704 Collaboration [45], at y/s = 19.4 GeV. These 
data cover a large range of (positive) values for xp, with px in the range [0.7 2.0] GcV/c. In what follows, we will 
use for our calculations an averaged value of pt = 1-5 GeV/c. Although this value might be slightly higher than the 
effective averaged pr of the data set, it is well indicative of the kinematic regime analysed. 

We start by noticing what the E704 data show: at large positive xp the SSA, A^, for charged pions is large and 
almost of the same magnitude, but opposite in sign for tt"*" and 7r~ (positive for 7r+). In order to describe these 
results in our approach invoking the Sivers effect alone we need a positive (negative) up (down) Sivers function, 
A'^ fq/pi {x, k±). This means that the non-leading flavour contributions (e.g. d t^^)) enter with opposite sign 
(all other terms in Eq. (44) giving positive contributions) with respect to the leading ones (e.g. u 7r+) in each 
charged pion SSA asymmetry. In principle this might reduce the magnitude of An and in some cases prevent a 
good description of data. This is particularly true for n~ production, where cancellations among different flavour 
contributions may be stronger. In this case, in fact, the dominance of the d — > 7r~ over the « — > 7r~ contribution 
in the fragmentation process can be weakened by the expected dominance (at large values of x) of the u over the d 
PDF's of the polarized proton. 

Let us also remind that the available sets of pion FF's are essentially for the neutral pion case, coming from fits 
to data with no flavour/charge separation. Therefore we will fix the parameters entering the Sivers function by 
reproducing at best the neutral pion SSA data; only in a second step, we look at Ajsi for charged pion production. 
This is slightly different from the procedure employed in former studies [7], where the used FF set were valence- like 
dominated (similar to our KKP-1 model set) and, adopting a simplified kinematics, a simultaneous best fit of charged 
and neutral pion SSA was performed. 

Within these assumptions we get for the three different FF sets discussed in Section II C 2, the following optimal 
sets of parameters for the Sivers function 

Nu = +0.40 a„ = 3.0 6„ = 0.6 

Afd = -1.00 ad = 3.0 6d = 0.5 (Kretzer) , (46) 



Nu = -1-0.40 a„ = 3.6 6„ = 0.6 

iVd = -0.55 ad = 3.0 6^ = 0.3 (KKP-1), (47) 

Nu = +0.40 a„ = 2.0 6„ = 0.3 

Nd = -0.90 ad = 2.0 ba = 0.2 (KKP - 2) . (48) 

In all three cases the parameter r entering the Gaussian k± dependence of the Sivers function is fixed at r = 0.7 (see 
also Ref. [20]). 

Note that we are employing the MRSTOl set [25] for the impolarizcd PDF's entering Eq. (43) instead of the GRV94 set 
[52], adopted for instance in Ref. [20]; as a consequence, the parameters for the Sivers function are slightly different. 
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Corresponding estimates for SSA are shown in Fig.s 15, 16, and 17. We first notice that using the Kretzer FF set, 
where the non-leading quark FF is relevant over a large range of z values, a simultaneous good description of neutral 
and charged pion SSA's seems to be more difficult (see Fig. 15). In the charged pion case, in fact, cancellations 
between leading and non- leading contributions in the numerator of Ajv are more effective. Even maximizing in size 
the Nd parameter, see Eq. (46), but keeping the agreement with tt" SSA data, we are not able to give a description of 
the Aat data as good as for the other two sets. Apparently, the best agreement with data is obtained with the KKP-2 
model set (see Fig. 17) where, without any need of maximizing the parameter, see Eq. (48), the concordance 
with neutral and charged pion SSA is quite good. Concerning the KKP-1 model set, where non-leading quark FF's 
are strongly suppressed, it also allows a reasonable description of data (see Fig. 16), being a factor of 2 far from 
saturating the positivity bound for both flavours, Eq. (47). Notice that in all cases the a, h parameters cannot be 
changed very much if one wants to correctly reproduce the shape of the SSA. Therefore, the above conclusions can 
hardly be changed by modifying these parameters. 

We have also considered a different functional form for Tl{k±) in the fc_L-dependent part of the Sivers function, 
Eq.s (39) and (42): 



n'ik,) = ^^, (49) 



2k±Mo 

where Mq is a typical hadronic scale parameter. While still ensuring the positivity bound, Eq. (40), this form is 
less (power-like) suppressed at large kj_ values as compared with the Gaussian-like form of Eq. (42) adopted in our 
calculations. We have explicitly checked that taking Mq = 1//3 = 0.8 GeV/c (a natural choice that also minimizes 
the number of free parameters), and keeping the same sets of parameters as given in Eq.s (46), (47), and (48), we are 
still able to describe the E704 data for An, at the same level of accuracy as shown in the plots discussed above. 

Another result of this analysis, which consistently treats the full k±_ kinematics in the calculation of both the 
numerator and the denominator of An, is that all the main features of the Sivers effect and Sivers function found in 
previous papers [7] based on a simplified kinematics, are confirmed. On the contrary, the inclusion of full partonic 
kinematics seems to substantially reduce the role of the Collins mechanism in the p^p tt X process (sec [19]) 

As a next step, we may now use these findings to give estimates of ^jv for pion production at the much larger 
energies reached at the RHIC collider, namely ^/s = 200 GeV, but almost in the same pr region (so that the so- 
far unknown pQCD evolution of the Sivers function is not relevant). In this sense the following estimates can be 
considered as real predictions of our approach. Indeed they were obtained in a preliminary version of this study [53] 
before the new STAR data for neutral pion SSA [54] became available. Note however that in this analysis we have 
not considered the possible role of Sudakov suppression factors [55]. 

In Fig. 18 we show our results obtained using the three sets of parameters for the Sivers function discussed above 
and compare them with STAR experinic^iital data [54]. Curves are calculated at fixed pseudo-rapidity rj = 3.8, so 
that pT and xf are correlated. To be safe with our pQCD approach, we consider only the region of pr >1-1.5 
GeV/c; our curves are then cut at pr values around 1.25 GeV/c (which implies xf ^ 0.3). This comparison indicates 
reasonable agreement with data, although error bars are still quite big at large xp. We stress the fact that, as 
predicted theoretically, sizable pion SSA, with features similar to those found in the E704 kinematical range, are 
confirmed experimentally even at such large cm. energies. 



B. Sivers effect in p'p — » jX 



Direct photon production in singly polarized pp collisions is definitely of interest in this context, since it allows 
(unambiguous) access to the Sivers effect. In this case there is clearly no contamination from fragmentation processes 
(Collins effect); in the positive, moderately large Xp region only the Sivers effect could be responsible for such an 
asymmetry in our approach. Therefore, the numerator of the SSA, Eq. (32), is given by an expression similar to that 
for the unpolarized cross section, Eq. (23), substituting the unpolarized PDF fa/p{xa,k±a) with the corresponding 
Sivers distribution fg^/pt {xa, k±a); the denominator is simply given by twice the unpolarized cross section. 

Unfortunately the few available data, at very small xf, are compatible with zero within the (large) errors. Therefore 
a comparison with our predictions is somehow inconclusive. No test of our assumption of a valence-quark dominated 
Sivers function is possible either. 

In Fig. 19 we show our predictions, utilizing the three different Sivers function parameterizations considered, for 
the kinematics of the E704 experiment, that is at \/s = 19.4 GeV and fixed pT = 2.7 GeV/c, as a function of xp- The 
calculated asymmetry is sizable at medium- large values of Xf, and could be measurable in principle. Experimental 
data at intermediate xp values would be of great help. 

We also show our results for STAR kinematics in Fig. 20, namely at ^/s = 200 GeV and fixed rapidity y = 3.8: 
once again sizable SSA are found which could be checked experimentally. 
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We can then conclude that SSA in direct photon production could be an important tool in understanding the Sivers 
mechanism. The same is true for any process where only one at a time of all possible effects might be at work. 
This is the case for instance of the single polarized DrcU-Yan process, where by suitably integrating on the angular 
variables identifying the leptonic scattering plane only the Sivers effect would survive [20] ; another key process could 
be heavy meson production in hadronic collisions, where the dominance of quark and gluon annihilation channels can 
enhance or even select out the Sivers effect [51]. 

IV. CONCLUSIONS 

In this paper we have studied the role of partonic intrinsic transverse momentum in inclusive particle production 
in (im)polarizcd high-energy hadronic reactions. Several papers in the framework of LO and NLO coUincar pQCD 
have already been devoted to this subject in recent years. Our analysis is different in that it is mainly aimed at 
stressing the connection with recent approaches to the problem of the large single spin asymmetries measured in 
inclusive pion production at different cm. energies at large positive xp and moderately large px- A lot of theoretical 
work has shown that LO pQCD with the inclusion of partonic intrinsic transverse momcntiim and a new class of 
leading twist, spin and fe_L-dependent PDF's and FF's is a promising approach to this problem. However, in those 
papers it was never shown in detail to what extent this approach is able to reproduce unpolarized cross sections for 
the same processes in similar kinematical conditions. Therefore, in the first part of this paper we have addressed this 
problem in a systematic way, with a full account of k± effects both in the PDF's/FF's and in the elementary partonic 
process. We have applied our approach to the Drell-Yan process and to inclusive pion and direct photon production 
in hadronic collisions in several kinematical situations of interest, performing a detailed comparison with available 
experimental data. Our results show that, with few noticeable exceptions, our LO approach (when complemented 
with proper NLO if-factors) is in reasonable agreement with a large set of experimental data for unpolarized cross 
sections. Therefore, this gives support to the validity of the same approach in the study of single spin asymmetries. 

In the second part of the paper, we have then applied our approach to the study of SSA in inclusive pion and 
direct photon production, with the inclusion of Sivers contribution alone (different contributions, like the Collins 
effect, in these or other processes have been or will be studied elsewhere). We have shown that the E704 experimental 
results on pion SSA can be reproduced with good accuracy by using physically reasonable parameterizations of the 
Sivers function for valence quark contributions (therefore, neglecting sea-quark and gluon Sivers functions in first 
approximation). 

These same parameterizations lead to predictions in agreement with measurements performed by the STAR Col- 
laboration at the RHIC-BNL accelerator at much higher cm. energies and similar xp and pt ranges. 

We have also shown that our analysis of SSA in terms of the Sivers effect and with a full treatment of k± effects 
is in good qualitative agreement and confirms all main results and conclusions of former studies performed keeping 
only leading contributions in k± and using a simplified partonic configuration. 

However, some additional comments and some words of caution, regarding our results on impolarizcd cross sections, 
are in order. The use of LO and NLO pQCD calculations (already in the coUinear configuration) in the study of 
relatively low pt data faces a number of formal problems and leads to some model dependent choices regarding e.g. 
the factorization scale and the strong coupling constant behaviour. Moreover, from the numerical point of view, the 
lower the pr values considered, the more relevant are these model dependences. 

Our results, which cover a range of pr values among 1-15 GeV/c for the cross sections and of a few GeV/c at 
most for the SSA, have been obtained by performing a number of choices that we believe very reasonable and are 
largely adopted in the literature. Admittedly, much work remains to be done, in particular on the more formal aspects 
of the approach. We can however conclude with good confidence that there is not presently any strong argument 
or evidence against the applicability of our approach to SSA calculations coming from unpolarized cross section 
results, in particular for the kinematical situations of relevance for SSA. Within the uncertainties inherent the pQCD 
approach discussed above (hopefully less relevant for the SSA which are expressed as ratios of cross sections) there 
is good agreement among our results and experimental data. Concerning those few cases where our approach clearly 
fails to reproduce the results, we want to notice the following: 1) We have considered a very large set of data in 
different kinematical situations and from both fixed target and collider experiments; 2) Even NLO pQCD in coUinear 
configuration faces several problems in fitting simultaneously the large collection of data considered here. 

Let us finally add some words concerning transverse hyperon polarization in unpolarized hadronic collisions. It 
was shown that our approach is in principle able to reproduce most of the striking features of the experimental data 
available; the role of the so-called polarizing fragmentation function was emphasized. However, a combined analysis 
of SSA and impolarizcd cross sections, like that performed here for the Drell-Yan process and inclusive pion, photon 
productions, is still lacking. This problem deserves a separate analysis which is under way and will be presented 
elsewhere. 
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APPENDIX A: FULL k± KINEMATICS 



We give here a detailed treatment of partonic kinematics with proper inclusion of transverse momentum effects, 
along the same lines of [56] . Let us consider the hadronic reaction AB CX in the A-B center of mass frame with A 
moving along the positive Z axis and fix the scattering plane as the X-Z plane. The 4-momenta of hadrons A, B, C 
read 



p'X = {Ea,0,0,P) p^s = {EbM-P) p'^ = {Ec,Pt,0,Pl) 



(Al) 



with Ea,b = y'P^ + Tn\^Q and Ec = \Jpt +pi + "^c- ^'^^ equal-mass initial hadrons [riiA — rriB = m), we have 

Ea = Eb = \/i/2, with s = (pa+Pb)'^- 

For massless partons a, b inside hadrons A, B we introduce light-cone momentum fractions Xa = p^ lv\ i = Pb /Pb 
and transverse momenta k±a, l^±b- Their four-momenta then read 



XaW— \ 1 



XbW- 



A-2 



2fc±a 



: cos ( 



XaWv'S 



sin (pa, 1 



xfiu'^s^ XbW^/s 



cos (pb, 7= sm (pb, —1 



XbW^ 



'^±b 
xlw'^s 



(A2) 



where the factor w, defined bsw = 1-F V1-4toVs /2, can be safely taken equal to one for high-energy processes, as 

considered in this paper. In Eq. (A2) fcxa,6 = |fe-La,6|) and (j)a,b are the azimuthal angles of parton a, h three-momenta 

in the hadronic cm. frame. 

The four- momentum of fragmenting parton c is given in terms of the observed hadron momentum p^, of the 
light-cone momentum fraction z = Pq/p^ and of the transverse momentum of hadron C with respect to parton c 
light-cone direction, k±c- In order to have all transverse momenta defined in the hadronic cm. frame, we define the 
(two-dimensional) transverse momentum, k±c, as a genuine three- momentum: 



kc = kc (sin 9k^ cos > sin 9kc sin cpkc , cos Okc ) : 



(A3) 



and impose the orthogonality condition kc ■ = 0, see below. Keeping hadron and parton masses into account, the 
parton four-momentum, = (EcPc), is then given as 



Pc 



E, 



m 



Vp'c - 



= {Pc - kc) = 



m 



c 



Vp'c - K 



= {Pt - kc sin Okc cos (pkc^-^c sin dkc sin (l)kc,PL - kc cos 0kc) A^^) 



Ec + Vpc - kl 



2z 



1 + 



{Ec + Vp^^) 



(A5) 



When, as it is the case in this paper, hadron C is a pion and only light quark and gluon fragmentation is considered, 
hadron and parton masses can safely be neglected in the kinematics; eqs. (A4) and (A5) then simplify to 



P. 



E, 



^ Ec + VEl-kl 

2z 



{Pt - kc sin 9kc cos (pkc , -kc sin sin (pkc ,Pl- kc cos 9kc ) 



(A6) 



(A7) 



On the other hand for heavy hadrons (like A hyperons) and/or heavy quarks (like charm quark fragmentation into D 
mesons), mass corrections can be comparable with fej_ effects. 
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Another ingredient entering our basic factorization formula, Eq. (4), is the Jacobian factor J{z,\kc\) connecting 
the parton c to hadron C invariant phase space, defined as 



-^ = -Jiz,\kc\)^, 



(A8) 



which for coUinear and massless particles reduces simply to J = 1. In the most general case, after some algebra, one 
gets 



J{z,\kc\) 



A{El-kl) 



(A9) 



where the last expression holds for massless partons and hadrons. 

The orthogonality condition kc ■ Pc = ^ mentioned above is realized through the following relation 

dPk^c = d^kc S{kc ■ pj , 



(AlO) 



where is the unit vector along the direction of motion of parton c. By using Eqs. (A3) and (A4), Eq. (AlO) becomes 

(All) 



d k±c = kcdkc dOkc d(j)kc 



\/Pc - fee 



PT sm v-^^,^ 



0^ Wkc - C) + 5{4>kc - (27r - C))] 



kc-PL COS Okc „ , , ^ 
^ PT sm 9ka --^kc- 



(A12) 



In this way the integration over (f)kc can be carried out directly (notice that there are two possible solutions to be 
considered) . 

With the expression of parton momenta given in Eqs. (A2) and (A6) one can calculate the partonic Mandelstam 
invariants; by exploiting the elastic, massless parton constraint s + £ + w = 0, one can in turn fix the value of z: 



T 

Z 

U 



S = {Pa+Pbf 

i = {Pa - Pcf 

U = iPb-Pcf 

s5{s + i+u) = zS{z + 



1 - cos((/)o - (j)b) 



A-2 i.2 



T + U 



(A13) 

(A14) 
(A15) 
(A16) 



where the two functions T and U, easily obtained from the explicit expressions of parton momenta, are independent 
of z. For heavy quarks, one can still fix z from the elastic constraint but in this case the expressions for t and u, in 
terms of z, are more involved. 

We now discuss, limiting to the massless case, the constraints on the full phase space entering our factorized 
expression for (un)polarized cross sections. Besides the trivial bounds < Xa,b, z < 1, < (j)a.b < 27r and < 6kc l£ t^, 
we require that, even including intrinsic transverse momentum effects, a) each parton keeps moving along the same 
direction as its parent hadron, Pa(b) ' b) > 0, and b) the parton energy is not larger than the parent hadron energy, 
Ea(^b) ^ -£'A(B)- This implies the following bounds 



k±a(b)/Vs < min Xa(b), y Xa(b){^ " Xa(b)) 



(A17) 



Analogously, for the fragmentation process c ^ C+X we require ■ Pc > and Ec < Ec (both fulfilled by Eq. (A7), 
where we have consistently disregarded the solution Ec = Ec — \J Eq — k^ /(2z)). The last constraint implies the 
following bound on kc, at fixed z 



kc/Ec <1 {z< 1/2); kc/Ec < 2^2(1 -z) (z > 1/2) . 
By requiring | cos^^^j < 1, see Eq. (A12), we have a further constraint on kc, at fixed 9kc, namely 

PL cos 0kc - PT sin 0kc < kc <Pl cos 9kc + Pt sin 0kc ■ 



(A18) 
(A19) 
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A word of caution is needed for the partonic flux factor when incoming partons are not coUinear. As shown in [57] 
the correct convolution formula connecting the partonic to the hadronic process is expressed in terms of partonic and 
hadronic tensors, that is in terms of scattering amplitudes squared rather than directly in terms of cross sections. 
The convolution then only involves parton-hadron light-cone energy ratios and no partonic flux, which would depend 
on the relative azimuthal angle between parton momenta, enters. 
More precisely, for the reaction AB — > cd one has 



/P+ p+ 
dXad'^k±adXbd^k±b fa/A{Xa, k±a)fb/B{Xb, k±b) \M\l^^^^ ^ ^ 
PaPb 

= (dXad^kji_adXbd^k±_bfa/A{Xa,kjLa)fb/B{Xb,k±^b)\M\lb^cd——- (A20) 

J XaXb 



In order to obtain the proper normalized AB cd cross section one has to divide Eq. (A20) by the hadronic flux (2s). 
On the other hand since our convolution formula is given in terms of partonic cross sections, which are normalized to 
the partonic flux 2s, to restore proper normalization we have to rescale them by the factor s/{xaXbs). 
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FIG. 1: Invariant differential cross section for the Drell-Yan process at ^/s ~ 19.4 GeV and fixed rapidity y — 0.4, as a function 
of the transverse momentum of the lepton pair qt and averaged over different invariant mass bins (see the legend). The 
parameterization MRSTOl [25] for the unpolarized parton distributions is used, with 1//3 = 0.8 GeV/c. Curves are rescaled by 
a fixed if -factor, K — 1.8. Data are from [27]. 
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FIG. 2: Invariant differential cross section for the Drell-Yan process at -^s ~ 23.8 GeV and fixed rapidity y = 0.21, as a 
function of the transverse momentum of the lepton pair and averaged over different invariant mass bins (see the legend). 
The parameterization MRSTOl [25] for the unpolarized parton distributions is used, with 1//3 = 0.8 GeV /c. Curves are rescaled 
by a fixed /S'-factor, K — 1.8. Data are from [27]. 
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FIG. 3: Invariant differential cross section for the Drell-Yan process at ^/s ~ 27.4 GeV and fixed rapidity y = 0.03, as a 
function of the transverse momentum of the lepton pair and averaged over different invariant mass bins (see the legend). 
The parameterization MRSTOl [25] for the unpolarized parton distributions is used, with 1/(3 = 0.9 GeV /c. Curves are rescaled 
by a fixed if-factor, K — 1.6. Data are from [27]. 
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FIG. 4: Invariant differential cross section for the Drell-Yan process at -^s ~ 38.8 GeV and fixed xf = 0.1, as a function 
of the transverse momentum of the lepton pair qr and averaged over different invariant mass bins (see the legend). The 
parameterization MRSTOl [25] for the unpolarized parton distributions is used, with 1/P — 0.95 GeV/c. Curves are rescaled 
by a fixed /S'-factor, K — 1.6. Data are from [28]. 
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FIG. 5: Diflterential cross section for the Drell-Yan process at ^/s = 62 GeV, as a function of the transverse momentum 
of the lepton pair gr and averaged over the invariant mass bin 5 GeV < M < 8 GeV and over the Feynman variable bin 
—0.1 < xf < 0.8, {xf = 1qLl\fs)- The parameterization MRSTOl [25] for the unpolarized parton distributions is used, with 
1//3 = 1.0 GeV/c. The theoretical curve is rescaled by a fixed isT-factor, K = 1.5. Data are from [29]. 
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FIG. 6: Invariant differential cross section for prompt photon production inpp collisions, at -^s ~ 19.4 GeV and fixed a;F = 0, as 
a function of the photon transverse momentum pr. The parameterization MRSTOl [25] for the unpolarized parton distributions 
is used, with 1//3 = 0.8 GeV/c (thick line). For comparison, the result in collinear partonic configuration (thin line) is also 
shown. Both curves are rescaled by a fixed if-factor, K = 1.6. Data [35] are averaged over the xf bin —0.15 < xp < 0.15. 
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FIG. 7: Invariant differential cross section for prompt photon production in pp collisions, at -^s ~ 23 GeV and different xp 
values, as a function of the photon transverse momentum px- The parameterization MRSTOl [25] for the unpolarized parton 
distributions is used, with 1//3 — 0.8 GeV /c (thick lines). For comparison, the results in collinear partonic configuration (thin 
lines) are also shown. All curves are rescaled by a fixed Tf-factor, K = 1.3. Data are from [36]. 
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FIG. 8: Invariant differential cross section for prompt photon production in pp collisions, at y/s — 63 GeV and fixed rapidity 
y = 0, as a function of the photon transverse momentum px- The parameterization MRSTOl [25] for the unpolarized parton 
distributions is used, with 1/P — 1.0 GeV/c (thick line). For comparison, the result in coUinear partonic configuration (thin 
line) is also shown. Both curves are rescaled by a fixed /T-factor, K = 1.5. Data [37] are averaged over the rapidity bin 
-0.2 <y < 0.2. 
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FIG. 9: Invariant differential cross section for inclusive tt" production in pp collisions, at different cm. energies (see the legend) 
and fixed rapidity y — 0, as a function of the pion transverse momentum pr- We use the parameterization MRSTOl [25] for the 
unpolarized parton distributions, with 1//3 = 0.8 GeV/c, and the parameterization KKP [40] for the unpolarized fragmentation 
functions, with /3' given in Eq. (24). Tf- factors for each case are shown (some data sets and their corresponding curves are 
further rescaled for clarity). Data are from [38, 41-44]. 
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FIG. 10: Invariant differential cross section for inclusive tt" production in pp collisions, at ^/s = 23.3 GeV and for two 
cm. scattering angles, as a function of xp. We use the parameterization MRSTOl [25] for the unpolarized parton distributions, 
with 1//3 = 0.8 GeV/c, and the parameterization KKP [40] for the unpolarized fragmentation functions, with (3' given in 
Eq. (24) (thick lines). For comparison, the results in coUinear partonic configuration (thin lines) are also shown. All curves are 
rescaled by a fixed Jf-factor, K — 2.5. Data are from [41]. 
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FIG. 11: Invariant differential cross section for inclusive vr" production in pp collisions, at ^/s = 52.8 GeV and for different 
cm. scattering angles (see the legend), as a function of xp- We use the parameterization MRSTOl [25] for the unpolarized 
parton distributions, with 1//3 = 0.8 GeV/c, and the parameterization KKP [40] for the unpolarized fragmentation functions, 
with /3' given in Eq. (24) (thick lines). For comparison, the results in collinear partonic configuration (thin lines) are also 
shown. All curves are rescaled by a fixed Jf-factor, K — 2 (some data sets with their corresponding curves are further rescaled 
for clarity). Data are from [41]. 
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FIG. 12; Invariant differential cross section for inclusive tt" production in pp collisions, at ^/s — 19.4 GeV and for different pr 
values (in GeV/c, see the legend), as a function of xf- We use the parameterization MRSTOl [25] for the unpolarized parton 
distributions, with 1//3 = 0.8 GeV/c, and the parameterization KKP [40] for the unpolarized fragmentation functions, with /?' 
given in Eq. (24) (plot (a), thick lines). For comparison, the corresponding results in collinear partonic configuration (plot (b), 
thin lines) are also shown. All curves are rescaled by a fixed AT- factor, K = 2.5. Data are from [42]. 
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FIG. 13: Invariant differential cross section for inclusive tt" production in pp collisions, at ^/s = 200 GeV and for two values 
of the pseudo-rapidity, 77, as a function of the pion energy E-^. We use the parameterization MRSTOl [25] for the unpolarized 
parton distributions, with 1//3 = 0.8 GeV/c, and the parameterization KKP [40] for the unpolarized fragmentation functions, 
with 13' given in Eq. (24) (thick lines). For comparison, the results in collinear partonic configuration (thin lines) are also 
shown. A"- factors for the two cases are shown. Data are from [46]. 
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FIG. 14: Ratios of invariant differential cross sections for inclusive charged pion production in pp collisions, at different LAB 
energies and fixed scattering angle 9 — 77 (mrad), as a function of pr- We use the parameterization MRSTOl [25] for the 
unpolarized parton distributions, with 1//3 — 0.8 GeV/c, and three different sets for charged pion FF's with the corresponding 
l3' parameters as in Eq. (24) and (25) (see text). Data are from [48]. 
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FIG. 15: An for inclusive pion production in pp collisions, at ^/s — 19.4 GeV and fixed pr ~ 1.5 GeV/c, as a function of xf- 
The parameterization MRSTOl [25] for the unpolarized parton distributions is used; fragmentation function set is from [47]. 
For the Sivers function, see Eq.s (41) and (43), parameters are given in Eq. (46), with 1//9 = 0.8 GeV/c and r = 0.7. Data are 
from [45]. 
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FIG. 16: An for inclusive pion production in pp collisions, at y/s = 19.4 GeV and fixed px = 1.5 GeV/c, as a function of xf- 
The parameterization MRSTOl [25] for the unpolarized parton distributions is used; fragmentation function set is KKP-1 (see 
Section II C 2). For the Sivers function, sec Eq.s (41) and (43), parameters are given in Eq. (47), with 1//3 = 0.8 GeV/c and 
r = 0.7. Data are from [45]. 




FIG. 17: An for inclusive pion production in pp collisions, at ^ = 19.4 GeV and fixed pr ~ 1.5 GeV/c, as a function of xf- 
The parameterization MRSTOl [25] for the unpolarized parton distributions is used; fragmentation function set is KKP-2 (see 
Section II C 2). For the Sivers function, see Eq.s (41) and (43), parameters are given in Eq. (48), with 1//3 = 0.8 GeV/c and 
r = 0.7. Data are from [45]. 
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FIG. 18: An for inclusive 7r° production in pp collisions, at ^/s = 200 GeV and fixed pseudo-rapidity rj = 3.8, as a function of 
xf- The parameterization MRSTOl [25] for the unpolarized parton distributions is used. Curves are for different fragmentation 
function sets and corresponding Sivers function parameterizations (see text). Data are from [54]. 
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FIG. 19: An for inclusive photon production in pp collisions, at ^/s — 19.4 GeV and fixed pr = 2.7 GeV/c, as a function of 
Xf- The parameterization MRSTOl [25] for the unpolarized parton distributions is used. Curves correspond to different Sivers 
function parameterization sets (see text). 
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FIG. 20: ^jv for inclusive photon production in pp collisions, at \fs = 200 GeV and fixed rapidity y = 3.8, as a function of 
xf- The parameterization MRSTOl [25] for the unpolarized parton distributions is used. Curves correspond to different Sivers 
function parameterization sets (see text). 



